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It is argued that the derivative expansion is a suitable method to deal with finite temperature 
field theory, if it is restricted to spatial derivatives only. Using this method, a simple and direct 
calculation is presented for the radiatively induced Chern- Simons-like piece of the effective action of 
(2+l)-dimensional fermions at finite temperature coupled to external gauge fields. The gauge fields 
are not assumed to be subjected to special constraints, and in particular, they are not required to 
be stationary nor Abelian. 
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I. INTRODUCTION 

Finite temperature field theory jjj is notoriously more 
difficult to deal with than its zero temperature counter- 
part. Lorentz symmetry is reduced and time and space 
play different roles. This is apparent in the imaginary 
time formalism, in which the time is Wick-rotated and 
compactified to a circle. As a consequence, tools such 
as the perturbative expansion or the derivative expan- 
sion, that were quite useful at zero temperature may 
become unreliable. These remarks apply, in particular, 
to (2 + l)-dimensional QED. In this theory, due to the 
compactification, topologically non-trivial (large) gauge 
transformations are supported at finite temperature [|2|. 
In principle, one would expect that the coefficient of the 
radiatively induced Chern-Simons (CS) term in the ef- 
fective action || would be correctly quantized Q, con- 
sistently with gauge invariance of the partition function 
at any temperature || . This turns out not to be the case 
when that coefficient is computed using perturbation the- 
ory (|, 0. Existing (Abelian and non- Abelian) exact re- 
sults for particular configurations [j], ||, [| show that 
gauge invariance is an exact symmetry but perturbation 
theory fails to see this; perturbation theory is designed 
to describe the effective action in the neighborhood of 
a vanishing gauge field and large gauge transformations 
necessarily move the gauge field configuration away from 
the perturbative region. The derivative expansion is also 
problematic as already shown in a (0 + l)-dimensional 
setting |lC| (using the real time formalism). The prob- 
lem is clear in the imaginary time formalism since the 
energy becomes a discrete variable. This suggests that 
the trouble exists only if one insists in taking a deriva- 
tive expansion in both time and space. Technically, to 
make a derivative expansion of the effective action func- 
tional W, means to replace the given configuration A^(x) 
by a family of configurations A^ (Xx) , where A is a book- 
keeping parameter, and then expand W(X) in powers of 
A plf. At finite temperature, this procedure is inconsis- 
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tent, since the external (bosonic) fields are required to be 
periodic in the (Euclidean) time and the dilatation in the 
temporal direction breaks this constraint. On the other 
hand, no problem should appear if the derivative expan- 
sion is restricted to the spatial directions (or at least, 
no new pathologies as compared to the zero temperature 
case). Moreover, gauge invariance can be accounted for 
if the bookkeeping parameter is introduced in the form 
Aq(x q , Xx), XAi(x , Xx). 

The exact results in |], 0, || were obtained for par- 
ticular configurations of the gauge field using specially 
adapted methods. In this Letter we show that general 
configurations are also amenable to explicit calculation 
using the method of expanding the effective action in the 
number of spatial covariant derivatives. We present a 
direct calculation of the induced CS term at finite tem- 
perature which yields a simple explicit expression for this 
quantity. This expression holds for general configurations 
and agrees with all previously known exact results. Of 
course, in addition to this, the effective action contains 
further contributions which are completely regular (i.e. 
gauge and parity invariant and ultraviolet finite) which 
are also calculable within the same general scheme. 

The fermionic effective action is defined as W — 
— logDetD, where D = + m is the Dirac operator, 

m is the fcrmion mass and = + A^ is the covari- 
ant derivative p2| . The Dirac operator acts on a single- 
particle Hilbert space which contains space-time, Dirac 
and internal degrees of freedom (i.e. those associated 
to the gauge group). The finite temperature T = 1//3 
is introduced by compactifying the Euclidean time. As 
usual the fermionic wave functions are antiperiodic and 
the gauge fields are periodic in the time variable. 

As is well known, the effective action is afflicted by sev- 
eral related pathologies, namely, ultraviolet divergences, 
many-valuation and anomalies in some of the classical 
symmetries of the Dirac operator. Our purpose is to iso- 
late precisely those pieces of the effective action which 
can have an anomalous or many-valued contribution. 
This excludes ultraviolet finite pieces, which are always 
regular. The expansion in the number of spatial co- 
variant derivatives (TTJ, [t2| is appropriate for this kind 
of calculations, since this expansion preserves gauge in- 
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variance order by order and, in addition, terms beyond 
second order are ultraviolet finite. A further simplifica- 
tion can be achieved by selecting only those terms which 
have abnormal parity, i.e., those containing a Levi-Civita 
pseudo-tensor, since the normal parity component of the 
effective action can be renormalized preserving all clas- 
sical symmetries and is one- valued |23| |, In passing we 
note that "parity" (defined as space reflection but in- 
cluding m — > — m) is a symmetry of the classical action. 
Therefore, we would expect the normal (abnormal) par- 
ity component of the effective action to be an even (odd) 
function of the mass. This symmetry may be anoma- 
lously broken as a consequence of gauge invariance, as in 
the massless case j|. 



II. THE CURRENT 

Following Schwinger p3| , we will work with the cur- 
rent, which is better behaved than the effective ac- 
tion. The current is defined through the relation SW = 
J d? x tr (J M 5 A (where SA^ is a local variation, J^x) is 
a matrix in internal space and tr refers to this space), 
and thus 
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(1) 



(tr/) refers to Dirac space.) In view of our previous 
remarks, our approach will be to compute the relevant 
terms of the spatial current Ji and then reconstruct the 
effective action from it. This requires to expand the cur- 
rent in powers of D; retaining only terms with and 
with precisely one Di. A suitable way to deal with the 
matrix element at coincident points, which combines well 
with gradient expansion, is to use the method of symbols 
Jlij , adapted to finite temperature [ pT[ [l^ ] and improved 
by Pletnev and Banin [ [l5[ [l6| . This gives 
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In this formula |0) is the constant wave function, i.e. 
(x|0) = 1 (which is periodic rather than antiperiodic) , 
Po = 27r(?i, + runs over the fermionic Matsubara 

frequencies, and finally, D M is related to the original co- 
variant derivative through a double similarity transfor- 
mation, 



D n 



dpi 



The inner similarity transformation yields D, t — > D„ 



(3) 



and corresponds to the original method of symbols. In 
that method, the integration over p cancels all contribu- 
tions except those where Di appears inside a commuta- 
tor (more precisely, in the form [Di, ]). In the improved 
method, this cancellation is achieved prior to momentum 



integration by means of the outer similarity transforma- 
tion in (ph. Explicitly: 
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where F^ — [D^D,,] and E{ = Foi is the electric 
field. The improvement does not extend to the com- 
pactified time coordinate, and indeed at finite tem- 
perature Do appears in two different ways which pre- 
serve gauge invariance, i) inside commutators, [Do, ], 
and ii) through the (untraced) Polyakov loop f2(x) = 



Texp(-j:° +(i A (t,x)dt) @. 



For our present purposes it will be sufficient to retain 
in (^) only those terms with at most one spatial covariant 
derivative. Upon substitution of these expressions in (0), 
and taking the Dirac trace keeping only the abnormal 
parity terms, produces 



J i (x) = 2imei 



-y 



d 2 p , 1 1, , 



O(Df). 



( 5 ) 

where A = m 2 + p 2 + (po — iD Q ) 2 . It is noteworthy that 
the result is directly ultraviolet finite and also consistent, 
i.e. a true variation p4| . 

The current (^|) picks up contributions from all terms 
with two spatial indices in the effective action, however, 
most of those terms are perfectly regular, in the sense 
that they are ultraviolet finite, one- valued, strictly gauge 
invariant and anomaly free, therefore we can simplify 
the calculation by isolating only the anomalous terms, 
namely 
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Jf n = imeij — 



(5^ 

po 



d 2 p 



(0) 



This "anomalous" current is still consistent. The differ- 
ence between the two currents, J~ — Jf n , involves [Do, Ei] 
and derives from (fully calculable) regular terms in the 
effective action which are quadratic in E- L |l7). On the 
contrary the anomalous current does not contain such 
a commutator and so it cannot derive from one of those 
terms. After carrying out the momentum integration and 
the sum over frequencies, the anomalous current takes 
the simple explicit form 



Jr = 6 u (x|{^(D ),D J }|0), 

where we have introduced the function 
i sinh(/3m) 



8-7T cosh(/3m) + cosh(/3a) 



(7) 



(8) 



This current preserves parity (as defined above), and can 
be written in a manifestly gauge invariant form by using 
the identity e'^ " = fl(x) p): 



sinh(/3m) 



E, 



2cosh(/3m) + n + n~ 17 3 



(9) 
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Note that this current gives the static limit of the 
gauge field self-energy instead of the long-wave limit 
JL8| . The latter limit follows from restoring the regular 
terms stripped from the anomalous piece. Our criterion 
for selecting the anomalous piece is that the remainder 
is strictly gauge invariant, however, regarding the self- 
energy after rotation to real time, the anomalous piece 
is analytical at zero momentum and energy and the non- 
analytical behavior is contained in the remainder jlTj]. 

At this step we can already verify that the zero tem- 
perature limit comes out correctly. Indeed, for large (3 



J i n = ^<™)eiJ E J > ( r = 0) 



(10) 



(where e(m) stands for the sign of m). This current de- 
rives from 



W an = -^e(m)W C s , (T = 0), 
where TUcs is the Chern-Simons action 



(11) 



U < s = ^ / d 3 xe tJ tv(A F tJ - Ad a A 3 ) . (12) 



Eq. ( p"T| ) is the standard induced Chern-Simons term at 
zero temperature ^ [| |l9) . 



III. THE EFFECTIVE ACTION 

Using the techniques developed in |ll], |l2| it is possible 
to reconstruct the anomalous piece of the effective action 
from its current while preserving manifest gauge invari- 
ance at every step of the calculation. This method will be 
presented elsewhere |l7j . A more direct approach, which 
we will follow here, is to proceed by fixing the gauge [p5| . 
We choose the gauge so that Aq(x) is stationary 



d A (x) = 0. 



(13) 



Such a gauge always exists (globally) Q. In this gauge, 
(||) becomes 



Jr(x) = e ij {tp(A (x)),E j (x)} 



(14) 



The final step is to recover the effective action associated 
to this current. To this end we split the current into two 
parts + J- 2 \ corresponding to the decomposition of 
the electric field into Ei = [Dq, Ai] — diA$. Noting that 
in our gauge A and D commute, it is readily verified 
that the first piece derives from 



WW 



d 3 xe i:j tr (ip(A ) A,[D , A 3 



(15) 



Thus the remainder must satisfy 

SW^ = - f d\n (SAt MA ), djA }) (16) 



To proceed, we further restrict the gauge (by means of a 
suitable subsequent stationary gauge transformation) so 
that Aq(x) is everywhere a diagonal matrix, i.e 

[Mx),diM*)]=0- (17) 
For subsequent reference we note that in this gauge 

[D ,[Do,Ai\} = [D ,Ei\. (18) 
With this choice of gauge, ( |l6| ) has the simple solution 

= 



d i xe ij tT{^>{A )d l A j ) 



provided that $'(a) = 2ip(a) or, in our case, 



(19) 



$(a) = — - tanh" 1 (tanh(/3m/2) tanh(/3o/2)) . (20) 
2tt(3 

Using the convenient notation of differential forms (for 
the spatial indices only), the final result obtained by 
adding up ( |l5| ) and (|l^), takes the form 



Wan = 



j dx tr ($(A )(A 2 — B) + V {A Q ) A[D , A] 



(where A = Aidxi and B = dA + A 2 = ^Fijdxidxj is the 
magnetic field). This action reproduces the spatial com- 
ponents of the current. Since any further contributions to 
the effective action cannot contain Ai, they would be of 
the form J dxo tr (/ (Ao)dAodAo) and this vanishes iden- 
tically in our gauge. 

A subtle, but important, point is that, depending on 
the space manifold and the gauge group, a given con- 
figuration Aq(x) may not be globally diagonalizable (i.e. 
with continuity at all points). Nevertheless, the "diago- 
nal" gauge ( |l7| ) always exists locally, that is, for each of 
the patches covering the space manifold. Without loss of 
generality we can assume that the support of the local 
variation in (jl^) takes place inside one of the patches, 
thus integration by parts in x is allowed and the cor- 
rect current is obtained from (|l9|). Because the current 
(|l4| ) is patch independent, this suggests that the inte- 
grand in (^l|) is also patch independent. As shown in 
the next section, this is actually correct, confirming that 
the functional W an is well defined. For simplicity, we will 
momentarily assume Aq(x) to be diagonal in a global 
gauge. 

Eqs. ( |l0| ) and (11) show that the calculation is con- 
sistent with the known zero temperature limit, and this 
can also be verified directly from IU an . Next we can study 
the behavior of this functional under gauge transforma- 
tions. We have to distinguish between the "gauge fixing" 
transformation needed to bring the original configura- 
tion to the gauge where A is stationary and diagonal, 
on the one hand, and the allowed gauge transformations 
which preserve the gauge conditions (|l^) and (|l^), on 
the other. Of course, gauge invariance of the partition 
function under the gauge fixing transformation cannot 
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be studied using our functional. Such gauge invariance 
follows from general arguments, e.g. by using the £- 
function version of the effective action Q. The gauge 
conditions and ([Tt]) are preserved by two kinds of 
gauge transformations, namely, stationary and discrete 
transformations [[llj . The first class is that of transforma- 
tions which are stationary and diagonal, and it is easily 
verified that W an is invariant under such transformations. 
The second class refers to those transforniations of the 
form U = exp(xoA), where A is a constant diagonal ma- 
trix with entries Xj = 2irinj/f3 (with rij integer). Under 
discrete transformations Ai transforms covariantly (i.e., 
homogeneously) whereas Aq — > Ao + A. The current Jl4| ) 
transforms covariantly (due to periodicity of <p(a)) thus 
Wan can only change by constant additive terms, that 
is, terms that do not give a contribution to the trans- 
formed current. This can be verified directly from the 
functional as follows. First, at a formal level the function 
<J>(a) is periodic and W an is invariant. More correctly, on 
its Riemann surface $(A + A) = $(A ) + ze(m)A/47r 
and this induces a variation —ie(m) /4ir J dxg tr(AcL4) on 
Wa n . This variation is precisely the same as that com- 
puted for the zero temperature functional in ( |i"l| ) un- 
der the same gauge transformation. This implies that 
this variation is of the form iiirkn, where k is an inte- 
ger depending on the gauge group and the gauge field 
configuration and n is the winding number of the gauge 
transformation |j26f . It can be noted that a perturba- 
tive expansion in W an or Jf n would break periodicity of 
the functions <p and $ and, in general, this spoils gauge 
invariance under discrete transformations. 



If there are large transformations and k is odd, the 
partition function changes by a factor (—1)" and gauge 
invariance must be enforced by adding a polynomial term 
i^Wcs @- Such term is present directly in the effective 
action when the ^-function regularization is used pi [l9| . 
This term is even as a function of m and so it introduces 
an anomalous breaking of parity. Note that this new term 
also adds a parity breaking contribution to the current; 
although the current (Q) is gauge and parity invariant it 
does not derive from a gauge invariant partition function. 
It also noteworthy that for massless fermions the factor 
(—1)" is also present || but comes from the normal parity 
component of the effective action, since W an vanishes in 
this case. 



I n H Hi it is shown that the full (abnormal parity) 
result is — J dxo tr (Q(Aq)B) provided that the configu- 
ration is stationary, Abelian and Aq is ^-independent. 
The same formula is found in Q in the non- Abelian sta- 
tionary case when the electric field vanishes. This is con- 
sistent with our W an noting that in both cases the Ai are 
stationary and commute everywhere with Aq, and so the 
explicit A cancel in (|2l|). 



IV. TOPOLOGICAL ISSUES 

Let us now analyze the general case in which the diag- 
onal gauge is not assumed to exist globally. In this case 
it will be convenient to work in the global gauge where 
Aq is stationary but not necessarily diagonal, and rewrite 
W an in that gauge. This yields 

W an = J dx Q tr (HA ){A 2 -B) + <p(A ) A[Dq, A]) 

(22) 

In this formula A denotes the vector field A of the di- 
agonal gauge rotated covariantly back to the stationary 
gauge. That is, if U denotes the stationary gauge trans- 
formation which brings Aq to a (local) diagonal gauge, 

A = A + dUU- 1 . (23) 

A crucial property of (^2|) is that it remains unchanged 
under the replacement 

A -> A' = A + C , (24) 

provided that the field C(x) satisfies 

d Q C= [Aq,C] =0. (25) 

This is readily verified noting that C is a 1-form and Aq, 
Dq and C are all commuting quantities. 

By construction A transforms covariantly (under 
transformations of the global stationary gauge and fixed 
local diagonal gauge), but is not globally defined since 
the gauge transformation U only exists locally. Let 
denote a set of local charts covering the space manifold. 
In each chart we can take a diagonalizing gauge transfor- 
mation and this defines a corresponding "covariant 
vector field" A^ . What have to be verified is that 2-form 
in (^2|) takes the same value in any of the charts. 

To show this let us note that the relation (|l8|), which 
was written in the diagonal gauge, in the stationary gauge 
becomes 

[Dq,[Dq,A^]] = [Dq,E]. (26) 

Therefore the field = A {k) - A (l) (which is station- 

ary due to (|H)) satisfies 

[A ,[A ,C^}}=0. (27) 

This alre ady implies that C^ k '^ commutes with Aq and 
satisfies (P5|). As a consequence any two A^ at the same 
point are related as in (^) and W an is patch independent. 

It is noteworthy that the field A' (x), denoting Aq in 
a diagonal gauge, is (or can be taken to be) equal in all 
patches. Indeed, the local eigenvalues of A (x), a\(x), 
can be labeled at each point so that they are continu- 
ous functions on the space manifold. Thus, without loss 
of generality we can choose the local diagonal gauges in 
such a way that the eigenvalues are ordered in all patches 
in the same way and so the A' (x), which are diagonal 
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matrices containing the eigenvalues, will be the same ma- 
trix in all patches. However, although A' (x) is global, it 
needs not correspond to a global gauge; in general, the 
eigenvectors corresponding to a\(x) cannot be chosen in 
a continuous way on the space manifold. Assuming the 
non degenerated case for simplicity, the choice of eigen- 
vectors in each patch will differ by a phase (and perhaps 
normalization, depending on the gauge group). These 
phases are contained in the diagonal transition matri- 
ces l/W U^ k \ which contain the non trivial topology of 
Aq{x) (under gauge transformations). 

When we bring a completely general gauge field con- 
figuration to an Ap-stationary gauge the result is not 
unique. The different possible stationary gauges so ob- 
tained are related either by a stationar y g auge transfor- 
mation or by a discrete transformation |11| . The discrete 
transformations are of the form U = exp(xoA), where 
A(x) is stationary, commutes with Aq(x) and has eigen- 
values Xj = 2iririj/l3 (with rij integer). The discussion 
of how Vy an changes under these transformations is es- 
sentially the same as that given before for the case of a 
global diagonal gauge. It can be noted, however, that the 
discrete transformations need not be topologically small 
if the diagonal gauge is not global (an explicit SU(2) ex- 



ample is shown in B, also discussed in JTT|.) 



V. CONCLUSIONS 

In summary, we have shown that the method of ex- 
panding in the number of spatial covariant derivatives 
is suitable to make computations at finite temperature, 
and in particular it is able to retain all subtleties usually 
tied to the effective action, such as topological pieces, 
many-valuation and anomalies. At the practical level, 
the method have been shown to be quite efficient, yield- 
ing a simple explicit form for the induced CS term at 
finite temperature, without assuming special constraints 
on the gauge field configuration. 
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